Spatio-temporal focussing of an ultrafast pulse through a multiply scattering medium 
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Pulses of light propagating through multiply scattering media undergo complex spatial and temporal distor- 
tions to form the familiar speckle pattern. There is much current interest in both the fundamental properties of 
speckles and the challenge of spatially and temporally refocussing behind scattering media. Here we report on 
the spatially and temporally resolved measurement of a speckle field produced by the propagation of an ultrafast 
optical pulse through a thick strongly scattering medium. By shaping the temporal profile of the pulse using a 
spectral phase filter, we demonstrate the spatially localized temporal recompression of the output speckle to the 
Fourier-limit duration, offering an optical analogue to time-reversal experiments in the acoustic regime. This 
approach shows that a multiply scattering medium can be put to profit for light manipulation at the femtosecond 
scale, and has a diverse range of potential applications that includes quantum control, biological imaging and 
photonics. 



INTRODUCTION 

The multiple scattering of coherent light is a problem of 
both fundamental and applied importance Q. In optics, phase 
conjugation allows spatial focussing and imaging through a 
multiply scattering medium |2-4|; however, temporal control 
is nonetheless elusive |5|, and multiple scattering remains a 
challenge for femtosecond science. When light propagates 
through a thick multiply scattering medium, the large number 
of scattering events exponentially deplete the ballistic photons 
with propagation distance, and the transmitted light is dis- 
persed in a highly complex manner. When illuminated by a 
single-frequency laser, such a medium is known to give rise to 
a spatial speckle | 9 |. Furthermore, this speckle is wavelength- 
dependent: illuminated by a broadband laser, the medium 
therefore also produces a spectral speckle whose character- 
istic size is inversely related to the confinement time in the 
medium, or Thouless time. It corresponds to a severe temporal 
broadening of the pulse |[TOl[TT|. Nonetheless, multiple scat- 
tering is deterministic and coherence is not destroyed. In the 
spatial domain, recent experiments have thus demonstrated 
that wavefront shaping is able to generate an intense focus 
through an opaque medium, and can be interpreted as phase 
conjugation |2-4|. For the acoustic and GHz-electromagnetic 
regimes, time-reversal experiments are the counterpoint to this 
principle in the time domain |6-8|: a short pulse incident on 
a medium generates a long coda, which may be time-reversed 
and retransmitted through the medium in order to regain the 
initial short pulse. These experiments have not only demon- 
strated that multiple scattering can be compensated but — 
more importantly — that this can actually improve address- 
ing, imaging or communication IE) [131. Due in particular 
to the inability to measure electric fields directly in the tem- 
poral domain at higher frequencies, an optical domain time- 
reversal experiment remains elusive; yet the ability to measure 



and shape femtosecond electric fields in the spectral domain 
nevertheless offers the opportunity of a route to the same goal. 

From a quantum-control perspective, the use of optimally 
shaped pulses to control light-matter interactions |14 | has 
been a fertile field of research 1 15 - 18|. The absence of pre- 
designed control mechanisms in complex systems may be ob- 
viated by a closed-loop feedback scheme 1 19 1 to find efficient 
pulse shapes at the expense of insight into the physical mech- 
anism |[T5l . By contrast, theoretically tractable model systems 
permit an open-loop approach where the appropriate control 
pulse is calculated directly |20, 21 1. In parallel, the propaga- 
tion of ultrafast pulses through strongly dispersive media has 
been a topic of research for many years. Atmospheric aberra- 
tions have been successfully compensated |22 | as well as the 
dispersion due to thick optical media 1231 . However, despite 
some insight into the multiple scattering of ultrafast pulses 
ll24l - [26l , measurement of the scattering medium transfer func- 
tion at a single speckle grain via scanning the frequency of a 
c.w. laser f27l[28J, and some two-photon excitation experi- 
ments in scattering samples |5 |, the effects of multiple scat- 
tering have ultimately been considered too complicated to be 
compensated. 

In this paper, we initially demonstrate a full spatio-temporal 
characterization of a femtosecond speckle field using spa- 
tially and spectrally resolved Fourier-transform interferome- 
try (SSI). Due to the linearity of the scattering process, knowl- 
edge of the spectral phase facilitates active temporal focussing 
of the speckle via the open-loop feedback of the measured 
phase to a spectral pulse shaper placed before the sample. 
Here we give the first experimental demonstration of this ef- 
fect. This experiment relies on a spatially resolved phase mea- 
surement, since the lack of large-scale spatial homogeneity in 
the speckle field prevents the spatial integration typical of con- 
trol experiments. An effective phase compensation is spatially 
localized according to the correlation length; thus a spatially 
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FIG. 1: The experimental set-up. The output of a laser oscillator 
(Osc) is divided by a beamsplitter (BS) to form an interferometer. 
One arm passes through a pulse shaper (PS) and is focussed onto the 
sample (S) by lens LI. The output speckle field is imaged by lens L2 
onto the spectrometer slit (2D-SM), which performs a measurement 
of the spectral intensity spatially resolved along the slit (upper inset). 
The reference arm is combined with an adjustable delay r and angle. 
The spectral pulse shaper is a folded 4/ line, and comprises a grating 
(G), cylindrical mirror (CM), plane mirror (M), folding mirror (FM) 
and spatial light modulator (SLM). 

resolved spectral shaper is not a prerequisite for a degree of 
spatial control. 



RESULTS 

Spatiotemporal characterization of the speckle field 

The schematic experimental layout is shown in Fig. 1 (see 
Methods for details). In essence, an ultrashort pulse passes 
through a pulse shaper and is focussed on an opaque sample. 
The sample is a scattering ZnO layer that transforms the ul- 
trashort pulse into a complex spatio-temporal speckle. This 
speckle is stationary on a timescale much longer than the ex- 
periment duration. The spatio-temporal speckle is reimaged 
and spatio- spectrally characterized in phase and amplitude via 
SSI (see Methods). 

Figure 2 shows the (a) intensity and (b) phase of a typ- 
ical spatio-spectral speckle field £^out(^,^) as measured by 
SSI. The spatio-spectral speckle is clearly demonstrated and 
the complex structure of E'out(^, y) is fully resolved in both 
phase and intensity — a prerequisite for quantum control. 



This structure is also visible in the one-dimensional spatial 
[(a), red, solid] and spectral [(a), green, solid] 'lineout' slices 
indicated. The integrated projections of the speckle inten- 
sity [(a), dotted lines], however, show that a spectrally unre- 
solved speckle image, as measured on a camera, would yield a 
strongly reduced contrast, while a non-imaging spectral mea- 
surement would only yield the initial source spectrum; this 
further motivates the necessity of a spatially resolved phase 
measurement for the temporal focussing experiment described 
below. Meanwhile the spectral phase reveals a similarly com- 
plex structure (Fig. 2b) ; as a consequence it is clear that a 
spatially averaged phase measurement offers no utility what- 
soever for pulse recompression. 

A Fourier transform of the complex field along the spec- 
tral axis gives the spatially resolved temporal behaviour of the 
speckle Eout{t,y) (see Fig. 3a). This spatio-temporal field 
exhibits the same complex speckle structure as before, as 
evinced by both the 3D plot and the spatial (red) and temporal 
(black) lineouts projected onto the walls. The spatially (black) 
and temporally (red) integrated fields are plotted above; the 
former reveals the confinement time which is fitted as approx- 
imately 2.5 ps, in good agreement with the spectral bandwidth 
measured from the autocorrelation function (Fig. 2c). For this 
function, the spatial and spectral correlation distances at the 
spectrometer were 50 /im (corresponding to 3.6 /im in the ob- 
ject plane) and 2.55 mrad/fs respectively. As a further remark, 
with an additional lateral scan, a full tridimensional speckle 
measurement could be envisioned. This would be interesting 
for fundamental studies of speckle properties such as vortices 
and singularities | 29 | . 

Focussing of the speckle field via spectral pulse shaping 

Multiple scattering is a complex yet linear process, there- 
fore the evolution of the pulse electric field through the 
medium may be elucidated by a Green function formal- 
ism. For a given spectral component, and for a distribu- 
tion of sources Es{r' ^ cj), the resulting field reads E{uj^ r) = 
JJJ Gscatt(^, r')Es{uJ, r')d^r\ where Gscatt(^, r, r') is the 
Fourier component of the Green function between point r and 
r' , in the presence of scattering. In our simple case where 
we enter with a focussed beam commensurate with the cor- 
relation distance of the scatterers Ei^{uj), the Green function 
reduces to Eout{^,y) = i^scatt(^, ^shaper (c^)^in(^) , where 

^scatt(^,^) = Ascatt(^,^)exp[z^scatt(^,^)] and i^shaper(^) = 

^shaper (^) Gxp ^shaper (^)] are the medium and shaper transfer 
functions respectively, with both consisting of phase and am- 
plitude contributions. For the measurements presented above, 
the SSI measurement reveals the relative phase 0scatt(^, y) — 
^in(^) with the pulse shaper set to i^shaper(^) = 1- The cou- 
pling between spatial and spectral modes in the scatterer, and 
the concomitant inseparability of i^scatt(^,^) into a product 
of spatial and spectral components, is the root of the exigence 
of a spatially resolved measurement for a spectral phase cor- 
rection. This coupling has very recently been exploited to ex- 
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FIG. 2: Experimental reconstruction of the spatio- spectral speckle electric field Eout{(jO,y) after propagation through a multiply scattering 
medium, (a) Speckle intensity \Eout{(^, y)\^, with respective projections onto the spatial (red) and spectral (green) axes — individual slices 
(solid lines) and integrated signals (dashed lines) are considered, (b) Speckle phase (j)oni{oj ^ y) , along with a spectral phase measurement 
0out(c(;, yo) localized along a single spatial slice, (c) Autocorrelation function | (£'out(^, y)^out(^^ y'))\^ demonstrating that the speckle is well 
resolved both spatially and spectrally. The spatial correlation length is related to the speckle grain size, while the spectral correlation length is 
the bandwidth of the medium, inversely proportional to the Thouless time. 



ert a degree of temporal control over isolated speckle grains 
via shaping of the spatial mode |[30l |3T1 . In general, a com- 
plete reversal of the scattering process would therefore require 
a spatio-spectrally resolved shaping element covering all the 
contributing modes, rather than successive shapers controlling 
the spatial and spectral degrees of freedom independently. 



In a subsequent experiment, such an open-loop compen- 
sation of the spectral phase is performed. A position is 
selected and the transfer function ^shaper(^) = — ^scatt(^, ^o) 
is applied to the pulse shaper. A successful phase compen- 
sation at yo is demonstrated by a further SSI measurement 
of 0out(^,yo) — ^inC*^, yo) = 0; furthermore, this flat output 
phase (and consequent temporal focus) is spatially localized 
as determined by the spatial phase correlation length. (It 
should be noted that phase-only shaping is sufficient to ob- 
tain a near-transform-limited pulse duration; amplitude shap- 
ing would serve predominantly to reduce the background in 
the temporal domain L32J ). 



DISCUSSION 

Experimentally, the shaped pulse is found to have a 
full width at half-maximum duration of 59 fs (close to the 
transform-limited duration of 54 fs); the optimal phase is 
reached after two iterations of the feedback procedure. In the 
temporal domain, this corresponds to an intense temporally 
focussed and spatially localized pulse that emerges from the 
background spatio-temporal speckle (Eig. 3b). The peak tem- 
poral intensity has a contrast ratio of 15 relative to the aver- 
age background before compensation, and the average back- 
ground along y^ is reduced by a factor of 2. As shown by 
the projected 'lineout' intensities, the temporal focus is spa- 
tially localized; the localization distance is 30 fim which is 
commensurate with the correlation distance of Eig. 2c. The 
spatially integrated temporal field (Eig. 3, top, black), exhibits 
the redistribution of the temporal intensity after shaping that is 
the signature of a temporal focus. The temporally integrated 
field, however, is not altered (Eig. 3 , top, red): the peak can 
be said to be spatially localized rather than focussed. Due to 
the spatially localized nature of the phase compensation, the 
peak integrated temporal intensity is not significantly altered. 
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FIG. 3: Reconstructed spatio-temporal intensities |^out(t, y)\^ (a) before and (b) after compensation of the phase at 2/0 = 2.66 mm. An intense 
peak emerges from the background with a contrast ratio of 15. One-dimensional 'lineouts' at the location of this peak (projections onto walls) 
and integrated signal (top) show that the peak is focussed in time (black) and localized in space (red). The temporal and spatial widths of the 
peak are the Fourier-limit pulse duration and the spatial phase correlation distance respectively. 



This ability to temporally recompress the pulse without spatial 
redistribution of energy may have important applications for 
the nonlinear imaging of biological samples, where a rise in 
fluence may result in the onset of thermally induced damage. 
As shown in time reversal and wavefront shaping experiments, 
another striking feature is that the more scattering the medium 
is, the more efficient the focussing will be. Indeed, the signal- 
to-noise ratio is governed by the ratio of the Thouless time of 
the medium over the initial duration of the pulse, i.e. the num- 
ber of independent spectral degrees of freedom. Furthermore, 
in contrast with conventional phase compensation techniques, 
here all shaping imperfections affect the signal-to-noise ra- 
tio; the temporal duration and spatial localization are limited 
by the Fourier-limit duration and autocorrelation size of the 
speckle respectively |[T3l . 

In conclusion, we have characterized the speckle spatio- 
temporal electric field of a multiply scattered ultrafast pulse. 
Furthermore, we have implemented an open-loop correction 
of the spectral phase in order to produce a naturally Fourier- 
and diffraction-limited pulse after the medium at a chosen po- 
sition. The results demonstrate the benefits of a spatially re- 
solved measurement: typical control experiments consider a 
spatial average, but for such systems that lack large-scale cor- 
relations such averaged measurements are valueless. More- 
over, we have demonstrated that these correlations permit spa- 
tial control without spatial resolution in the spectral pulse- 
shaper. Our study branches out previous spatial- speckle shap- 



ing techniques to the temporal domain and bridges the gap 
with time-reversal experiments in acoustics and electromag- 
netism. This capacity to recover a short pulse on a selected 
spatial speckle point has important potential applications for 
quantum control and photonics, as well as for the fundamental 
studies of complex media. It suggests that ultrafast diagnostic 
techniques including nonlinear microscopy and time-resolved 
spectroscopy may be performed deep within or beyond bio- 
logical tissue, beyond the range of ballistic photons. 



METHODS 
Experimental system 

The ultrafast source used for this experiment is an 80 MHz 
oscillator that delivers 4 nJ pulses at 800 nm with a spectral 
bandwidth of 80 nm. The oscillator pulse train is divided to 
form the two arms of an interferometer. In one arm is placed 
both a spectral pulse shaper and the scattering sample; the 
other arm acts as a reference and is recombined with an ad- 
justable relative angle and delay at the entrance slit of a home- 
built imaging spectrometer |[33]| . The spectral calibration was 
performed using a calibration lamp and verified against a com- 
mercial spectrometer; to ensure consistency the same calibra- 
tion was used for the pulse shaper. The spectrometer response 
— calibrated against a white-light source — varied by less 
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than 10% over the bandwidth of the pulse. The spectrome- 
ter has spatial and spectral optical resolutions of 20 /im and 
0.4 nm respectively; the two-dimensional detector permits the 
speckle spectrum to be spatially resolved along the slice that 
falls on the entrance slit (see inset of Fig. 1). The spatio- 
spectral intensity adhered to the well-known exponential de- 
cay probability distribution |9|, confirming that the speckle 
pattern was fully resolved. For a typical measurement 8 x 10^ 
oscillator pulses are integrated. 

The phase and amplitude pulse- shaper comprises a pair of 
liquid-crystal spatial light modulators (SLMs) in the Fourier 
plane of a folded double-pass 4/ zero-dispersion line 1341 . In 
order to optimize the trade-off between bandwidth and shap- 
ing resolution, the 4/-line grating is chosen so as to over- 
fill the SLMs; thus the output bears a clipped 30 nm hyper- 
gaussian spectrum and may be shaped arbitrarily within a 
23 ps time window. This output is then focussed to a waist 
of approximately 15 /im onto the surface of the sample. The 
scattered, transmitted light is collected by a lens with a nu- 
merical aperture of 0.25 and imaged onto the spectrometer slit 
with a magnification factor of 14. 



Sample preparation 

The scattering medium is thick layer of ZnO powder 
(a widely used white paint component) deposited homoge- 
neously on a microscope slide by sedimentation. The thick- 
ness has been measured to be L = 35 ± 5 /im, the transport 
mean free path has been measured asr = 2.1±0.2 /im, and 
absorption is known to be negligible. Since L > 10/*, the 
multiple scattering regime applies and virtually no ballistic 
light traverses the medium. 



SSI Technique 



The SSI technique fTSl |35l performs a relative measure- 
ment of the spectral phase between a reference and unknown 
pulse. The interference between the sample image and ref- 
erence beams causes interference fringes on the spectrometer 
whose spectral and spatial periods are determined by the rel- 
ative delay and angle of the two beams respectively. These 
fringes are additionally modulated by the relative phase be- 
tween the two beams. 



\A{uj,y)\'^\A{uj,y)\' 

\A,{uj,y) \ \ A,{uj,y)\ cos [0s -M^^v) - - 

(1) 



kyy] , 



Here r is the time delay between the two pulses and ky 
is the difference between the transverse components of the 
propagation vectors (such that their subtended angle is = 
ky / As, Ar, 0s and 0r denote the spatio- spectral ampli- 
tude and phase of the scattered (s) and reference (r) pulse re- 
spectively. 

In order to recover the amplitude and phase of the unknown 
pulse, the raw interferogram is thus Fourier transformed along 
both the spatial and spectral axes, and the relative spatio- 
spectral phase may be extracted through the filtering of an a.c. 
term followed by the inverse Fourier transform. This isolates 
one of the complex exponential terms that correspond to the 
cosine of the final summand of Eq. [TJ the argument of this 
exponential reveals the relative spectral phase modulo ±7r. 
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